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Interpretations of gamma-ray burst spectroscopy 
I. Analytical and numerical study of spectral lags 
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Abstract. We describe the strong spectral evolution that occurs during a gamma-ray burst (GRB) pulse and 
the means by which it can be analyzed. In particular, we discuss the change of the light curve as a function of 
energy and the spectral lag. Based on observed empirical correlations, an analytical model is constructed which 
is used to describe the pulse shape and quantize the spectral lags and their dependences on the spectral evolution 
parameters. Using this model, we find that the spectral lag depends mainly on the pulse-decay time-scale and that 
hard spectra (with large spectral power-law indices a) give the largest lags. Similarly, large initial peak-energies, 
Eq, lead to large lags, except in the case of very soft spectra. The hardness ratio is found to depend only weakly 
on a and the hardness-intensity-correlation index, 77. In particular, for low Eq, it is practically independent, and is 
determined mainly by Eq. The relation between the hardness ratio and the lags, for a certain Eq are described by 
power-laws, as a varies. These results are the consequences of the empirical description of the spectral evolution 
in pulses and can be used as a reference in analyses of observed pulses. We also discuss the expected signatures 
of a sample of hard spectral pulses (e.g. thermal or small pitch-angle synchrotron emission) versus soft spectral 
pulses (e.g. optically-thin synchrotron emission). Also the expected differences between a sample of low energetic 
bursts (such as X-ray flashes) and of high energetic bursts (classical bursts) are discussed. 
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1. Introduction 

The analysis of the prompt emission from gamma-ray 
bursts (GRBs) is giving us valuable clues to the envi- 
ronment from which the radiation emanates. Much evi- 
dence points to the fact that the gamma radiation stems 
from dissipat ion processes in a relativistically expan ding 
plasma wind ijRees fc Meszaroslll992HMeszaro£ll2002j) . ei- 
ther as shocks or in magnetic reconnections. Complex and 
variable light curves, that consist of several, often overlap- 
ping spikes or pulses, tell us through their variability and 
spect ral power distribution about the energetics ( K umaii 
19991), radial structure of the fireball ijBeloborodovl 
2000"), and can even give hints about the progenitor 
t,Kob avashi et al. (.200_2j; Salmonson & Galama (2002)). 
Individual pulses, which are assumed to be related to sin- 
gle emission episodes, give us information on the micro- 
physics; radiation processes (|Piran, 1990), and comoving 
properties of the fireball, such as th e densities and mag- 
netic fields ijRvde fc Petrosianll2002(l . We will here review 
the various spectroscopic analysis methods that are used 
to characterize the spectral evolution during bursts and es- 
pecially during their constituent pulses. We will examine 
the analysis of the spectral lags between the light curves, 
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measured in different energy bands, and how these results 
compare to the high-resolution spectroscopic analysis in 
which the time-evolution and correlations of spectral pa- 
rameters are deduced. The latter approach is useful as 
it can be used directly to test models of the dynamics 
and the emission processes. However, it can only be per- 
formed on cases in which we have high-resolution spec- 
tra and know their time-evolution. On the other hand, 
the lag description has the advantage that it can be mea- 
sured on practically every GRB, even weak ones, from 
all satellite missions which have more than one spectral 
channel. Moreover, several relations have been identified 
between the spectral lag and the physical quantities of 
the burst, such as the isotropically equivalent luminosity 
IINorris. Marani. fc Bonnell 1 (EOOO)) and the peak energy 
of the spectrum (|Amati et al. (2002)). It is therefore of 
interest to understand how these two descriptions relate 
to each other, since this will facilitate the interpretation of 
the spectral lag, and its dependence on other parameters. 

The observation that the light-curves in different en- 
ergy bands lag behind each other is a common feature 
in astrophysical objects, and is not only found in GRBs. 
There are several physical scenarios in which this can be 
explained. Positive lags, that is, the soft radiation lag- 
ging the hard, which is the dominant behavior observed 
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in GRBs, can be due to an intrinsic cooling of the radiating 
electrons, which will cause the radiation to dominate at 
lower and lower energies. Alternatively, a Compton reflec- 
tion of a medium at a sufhcient distance from the initial 
hard source will also cause a lag. Also, a convex surface 
that emits at relativistic speeds will cause the radiation 
emitted off axis, to be delayed and soft ened, thus produc- 
ing a lag, that is, the curvature effect ijRvde fc PetrosianI 
Furthermore, if there were an intrinsic time-scale 
that is constant from burst to burst, the change in viewing 
angle between the object and the observer, will cause both 
the observed time-scale and intensity to vary in a specific 
way due to t he an g ular dependence of the D oppler boost 
l)SalmonsonL l200d: lloka fc Nakamural l200l[l . A possible 
scenario for negative lags, that is, the hard radiation lag- 
ging behind the soft, is a hot medium, say a lepton cloud 
surrounding a cooler emitter, which will cause the soft ra- 
diation to be upscattered by inverse Comptonization. The 
photons get harder the more scatterings they suffer and 
thus they are more delayed. A definite answer to the rea- 
son for the spectral lags in GRBs has not yet been given. 

We will expand the spectral-lag analysis done sofar 
(den oted h e re as low-resolution spectroscopy, LRS: see, 
e.g., iBandl l)l997l) and references therein) by connect- 
ing it to the detailed spectral-evolution description that 
can be done for strong pulses (den oted here a s high - 
resolution spectroscopy, HRS: see, e.g., iFord et alJ lll995t): 
ICrider et all ^99^; iRvde fc SvenssonI ll2002t) ~he HRS 
and LRS are reviewed in Sect. El and in Sect. 01 we 
give an analytical treatment assuming a simple pulse 
model. The general empirical behavior of bursts, found 
through the observations made by the Burst and Transient 
Source Experiment (BATSE) on the Compton Gamma- 
Ray Observatory and their parameter distributions, are 
used in combination with the analytical findings to make 
a numerical study in Sect. 01 We simulate realistic pulses 
and determine how various spectral evolutions are mani- 
fested as spectral lags. Finally, we discuss the interpreta- 
tions that can be made in Sect. 5. 

In a subsequent paper II (Ryde et al. 2004), we will 
demonstrate the results on a sample of GRB pulses ob- 
served by BATSE and discuss the spectral lag correlations 
and various model scenarios. 



2. GRB Spectral Evolution: Empirical View 

2.1. High Resolution Spectroscopy (HRS) 

For the strongest BATSE pulses, spectroscopy is possible 
with high resolution, allowing us to deconvolve the ob- 
served count spectra through the detector response, thus 
providing the incoming photon spectrum. The deconvolu- 
tion is done through a forward-fitting method. This is done 
for several time intervals during the pulse, which allows 
details of the energy-flux spectrum to be followed in time 
and thus characterize the spectral e volution. Such stud- 
ies have been ma de for instance by IPreece et al.l l(200r)|) 
and iRvde fc Svensson (.2002) . However by necessity, the 



time-resolution becomes low and only the strongest and 
longest pulses can be studied in thi s way due to t he nee d 
for high signal-to-noise ratio, SNR I Preece et alJ l)l998l) V 
Using this method, the spectral evolution has been quan- 
tified in correlations between the observable parameters. 
The complete spectral and temporal evolution of a pulse 
can be characterized by three main observables, the peak 
of the EFe spectrum, Ep{t), the instantaneous energy 
flux, F(t), and the derived quantity, the energy fluence 
S{t) = J Fdt'. The relat ions between these three observ- 
ables are given by two empirical correlations: 

(i) the hardn e ss-int ensity correlation (HIC; see 
iBorgonovo fc Rvdel l(200ll) 'l. which relates the flux of the 
source to the 'hardness' of the spectrum (here represented 
hy E = Ep). For the decay phase of a pulse the most 
common behavior of the HIC is 



F = Fo{Ep/Eor, 



(1) 



where Eq and Fq are the initial values of the peak en- 
ergy and the energy flux at the beginning of the de- 
cay phase in each pulse and rj is the power law index. 
iBorgonovo fc Rvdd l)200j) studied a sample of 82 GRB 
pulse decays and found them to be consistent with a power 
law HIC in at least 57% of the cases and the power law 
index, r/ was found to have a broad distribution peaking 
approximately at 2.0 : rj = 2.0 ± 0.7. It should be noted 
here that even though the fits in some cases are not nec- 
essarily statistically consistent with a power law, in prin- 
ciple all cases do have a hard-to-soft trend that can be 
approximated by a power law. Indeed, by simply lower- 
ing the required lev el of s ignificance of the power-law fit, 
IBorgonovo fc RvdJ l)200l|) showed that even 75 % of the 
pulses were consistent with a power law. lRvde fc SvenssonI 
l|2002l) did derive alternative analytical descriptions of the 
HIC for pulses that do not follow exactly a power-law. 
Referring to their Fig. 1, it is however obvious that the 
power-law description can indeed be used as an approx- 
imate description. Below, we will therefore use the ana- 
lytical expression, given by Eq. ^ to describe the general 
behavior of all pulses. This description will then be exact 
for 60-70 % of pulses. 

(ii) the hardness-flu ence correlation, i?pk(f) (HFC; 
iLiang fc Kargatil l)l996(l ') which relates the hardness to 
the time running integral of the flux, that is, the fluence. 
Equivalently, the HFC represents the observation that the 
rate of change in the hardness is proportional to the lu- 
minosity of the radiating medium (or, equivalently, to the 
energy density); 

This corresponds to a linear relation between the Ep and 
the time integrated energy flux, the energy fluence. 



(3) 



where $0 is the decay constant. This behavior is most of- 
ten seen over the entire pulse and occurs in a vast major- 
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ity o f pulses ijLiang fc Kargatil Il996t iRyde fc SvenssoiJ . 
120021) . 

Ryde fc Svensson (2000) showed that these two cor- 
relations (Eqs. mni) form a complete description of the 
spectral evolution and define the light curve as well. In the 
original description they showed that the observed photon 
flux, F{E) /E or F(j/) /i/, can often be described by a recip- 
rocal function in time. In Sect.|21below we will reformulate 
their results in terms of the energy flux. 

2.2. Low Resolution Spectroscopy (LRS) 

To utilize the full temporal resolution of, for instance, the 
BATSE data, the deconvolution is not made and the in- 
coming light-curves are approximated by the count light- 
curves. This is also the case for weaker pulses for which the 
analysis described in the previous section is not possible. 
The spectral evolution is then instead quantifled through 
the change in shape of the pulse, in two or more broad en- 
ergy bands. For instance, the four channels of the BATSE 
discriminator rates (64 ms time resolution) can be used, 
and Beppo-SAX, 77ete-///FREGATE, and Swift have sim- 
ilar capabilities. Evolution of the spectrum, such as vari- 
ations in Ep and of the power-law slopes, will cause the 
light curves in the different energy-channels to be differ- 
ent. 

To quantify this change one can measure the shift in 
time between light curve pulses (spectral lags) and their 
change in width. These can be measured directly by using 
an analytical prescription of the pulses, but more often 
they are measured w i th th e use of correlation functions 
(e-g.- lFenimore et"al] l|l995l) ). The spectral lag. At, is de- 
fined as the lag, s, at which the cross correlation function 
(OFF) between the light curves in two different spectral 
channels, g and h, has its maximum: 

/oo 
g{T + s)h{T)dT. (4) 
'OO 

If the two light curves are very similar to each other, then 
the maximum occurs at a lag that measures the time that 
one is shifted compared to the other. For a detailed dis- 
cussion onthe general characteristics of CCF, we refer 
to lBandl lll997l) who gives the analytical expression of the 
CCF for a pure exponential decay (see his Fig. 1). Another 
way of quantifying the averaged, spectral evolution of a 
burst, using the LRS data, is to measure the hardness- 
ratio (HR) of the fluences (the time-integrated fluxes) in 
two channels. The HR and the CCF are not independent 

an d their relation is often soug ht. 

iNorris. Marani. fc Bonnelll l)200(]|) (see also iNorrij 
lj2002|) ) demonstrated that there is a trend that the HR is 
anti-correlated with spectral lag. A similar trend also ex- 
ists for the peak flux, in that the long-lag bursts all have 
low peak fluxes. Furthermore, the distribution of lag s is 
dominated by short lags, as also found by Band l)l997j) . A 
first step to understanding the connection between the de- 
tai led, sp ectral evolution and the spectral lags was taken 
bv lKocevski fc Liang. (,2003^ . They found that the HFC 



parameter $o has a nearly linear correlation with lag. An 
even tighter correlation was found when $0 was normal- 
ized to the peak flux. Fina lly, the HR and t he Ep are nat- 
urally correlated. Indeed. iBand et al.l l)2004() argued for an 
empirical relation where Ep is proportional to the square 

of the HR. 

IWu fc Fenimore I l)200(]l) pointed out that it is not al- 
ways reliable to determine lags with CCF. This is espe- 
cially the case for multi-peaked bursts, for which both the 
HR and the CCF give an average quantitative description. 
Variations of a burst's spectrum occur both within a pulse 
as well a s betw een pulses. This was conclusively shown 
bv iBandl l|l997l) . who used correlation functions and was 
able to conclude that ~ 90% of pulses have a hard-to-soft 
trend, while ^ 80% of the multi-peaked bursts have a soft- 
ening trend between peaks. Varying Eg, rj, and a between 
individual pulses will be important in determining the in- 
tegrated spectrum, which is the spectrum that is quan- 
tified by the lags and the HRs. \Wu fc Fenimore I l^200(t) 
also noted that the methods used to calculate the CCFs 
can affect the results considerably. For instance, the in- 
clusion of time intervals when the signal is at background 
will clearly affect the measurements. This adds complica- 
tions for the interpretation of these quantities. Below, we 
concentrate our study to an analytical model describing a 
single emission episode, resulting in an individual pulse. 
A comparison of the CCF measurements and the actual 
lags is further investigated in paper II. 



3. Analytical Model 

To catch the essence of the spectral evolution, we use the 
following simplified model: the emission episode consists 
of a single pulse which is assumed to be dominated by the 
decay phase, called a FRED pulse (fast-rise-exponential- 
decay). We also assume that the pulse has a monotonic 
decay of the spectral hardn ess, as measured b y Ep, that 
is, it is a hard-to-soft pulse ijFord et al.l l)l995(l ). The mo- 
tivation to study such a model is, first, that it is the indi- 
vidual pulses that bear the important physical signatures 
of the emission processes and these signatures are revealed 
mainly in the decay phase (^Rvde fc Petrosian. (,2002,) ). 
Most obs erved pulses do indee d ha ve such a sha,pe, a s 
noted by iKocevski etaD l)2003|) and iRvde et all l|2003l) . 
Second, most detailed spectral evolution studies have been 
made on the decay phase, since it often constitutes a major 
part of the pulse. Third, apart from the fact that many 
pulses are observed to be hard-to-soft pulses, pulses for 
which Ep is observed to track the fiux could very well 
still be intrinsicall y hard- t o-soft . A motivation for this was 
given by Kocevski et al.l l|2003l) : pulses that are emitted 
prior to the analyzed one often overlap it to some extent. 
Since, in general, the spectra of pulses soften with time, 
this contamination will contribute mainly soft photons. 
This will cause the measured spectrum to have an Ep 
that is lower than that of the spectrum of the analyzed 
pulse itself. Below, we derive the behavior of the decay 



4 



Ryde: GRB Spectroscopy 



phase of a pulse and make a few analytical estimates of 
the dependence of the lag on various variables. 

3.1. Decay Phase of a Pulse 

An analytical expression for the shape of a pulse can be de- 
rived based on the knowledge of the spe ctral evolution. We 
follow the calcula tions first outlin e d inlRvde fc Svenssoi] 
l)2000() (see also iKocevski et"all 1^200^ ) to derive the 
energy-flux decay of a pulse. The combination of the two 
empirical relations in Eqs. and Q gives the following 
differential equation governing the spectral evolution 

^0 ™ 



where Ep is the time derivative of the peak energy, Fq and 
Eq are the values of the flux and the peak energy at the 
start of the decay phase. The solution gives the energy 
flux 

F(t) = Fo l + ^V^ (6) 



T 



and the corresponding peak energy 
T 



Ep{t)^Eo 1 



where we have defined 
T = <^>oEq/Fo. 



(7) 



(8) 



To analyze the behavior of this function in the case rj = 
1, we note that the denominators above both have the 
appearance of 



1 



tjTV 
n J 



.with n = l/(l — ?]), 



(9) 



This expression has a limiting value of exp(— t/T) as n = 
(r/ — 1)^^ — > oo. Therefore as t] tends to 1, Eqs. ® and 
lO become 

F{t) = Foe-*/^ 
E,e~'/^. 



Epit) 



(10) 
(11) 



This shows that the full behavior is described by Eqs. 
© and O for all 77. We also note that it is possible to 
measure $0 directly from the light curve by identifying 
the time scale T — ^qEq/Fq, instead of by the analysis of 
the spectra otherwise necessary. 

Introducing d = rj/ (rj—l) (c? as in the asymptotic rfecay 
of the energy flux) we describe the peak and the energy 
flux decays as 



F{t) ^Foi^l- 
Ep{t) - Eo (1 



t 



T{d-l) 
t 



(12) 



(13) 



In Sect. I4.ll we discuss an expansion of these calculations 
to also include the rise phase, since we will need the full 
description in the numerical study. 



3.2. Spectral Lags 

Assume now that the two spectral channels between which 
the lag is measured, are characterized by the energies Ei 
and E2 . Assume further that the dynamical range of the 
change in Ep is larger or of the same order as Ei/E2. 
These assumption s are r easonable, which can be noted in 
iRvde fc SvenssonI ((20021 ^ . who studied the decay of many 
pulse decays and found that the typical dynamic range 
was ~ 4 — 5 (see their Fig. 1). Including a monotonic en- 
ergy decay during the rise phase of a pulse, this range 
will increase. Furthermore, typical values for BATSE are 
El ~ 25 keV and E2 = 100 keV. With these assumptions 
a portion of the non-power-law, curved photon spectrum 
will pass through both channels. The pulse shapes (peak 
and width), detected in the two channels, are then de- 
termined by how the spectral distribution changes with 
time; the main change of the spectrum is in the decay of 
Ep. The spectrum will then pass through a certain chan- 
nel at a speed dEp/dt which usually changes with energy. 
The spectral lag can be thought of as the time it takes 
for the spectrum to move from Ei to E2, which clearly 
is determined by dEp/dt. Note that the actual peak, _Ep, 
does not necessarily need to pass through the channels. 

Before discussing the analytical consequences of this, 
we will study three conceptually simple situations, which 
are valuable in illustrating the discussion. First, consider 
a situation where there is no spectral evolution, that is, 
the spectrum remains the same during the change in in- 
tensity-'^; F (X EIT^ with r] = 00. Then there will be no dif- 
ference in the light curve between the bands, apart from a 
normalization, and the spectral lag would be zero. Second, 
consider a situation where we allow for a spectral evolu- 
tion, but approximate the spectra with a Dirac (5-function 
in ene rgy (see, e.g., ICenI ()l999|) and iRvde & SvenssonI 
l)l999l) ) so that F-E.{E,t) = F(t) ■ d[E - Ep(t)], that is, 
with a monochromatic spectrum. In this situation, it is 
clear that the time interval between the pulse peaks in the 
two channels is exactly the time that the spectrum takes 
to move from Ei to £'2 . A third situation, that is still more 
realistic, is to take an energy spectrum that is described by 
a Heavyside function, Fe{E, t) = F{t) • e[£' - Ep{t)], that 
is, the spectrum consists of two power laws with photon 
indices a = —I and P = —00 that are sharply joined at 
Ep. As the break of the spectrum passes a certain energy, 
the flux disappears, which can constitute a characteristic 
point in the light curve. The difference in channel light- 
curves again will clearly reflect the Ep evolution. 

Based on these discussions, we therefore make the fol- 
lowing simplifled model that captures the general behav- 
ior, namely that the lag. At, is determined by 



dt 
AE dEp 



dE 



(14) 



and the rate, dEp/dt, is determined by the spectral evo- 
lution parameters according to Eq.((7|. Note also that, ac- 
cording to Eq. jSJ this rate is proportional to F/$o (see 

^ This situation never occurs for observed GRB pulses 
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also iKocevski fc Lian To describe this analyt- 

ically, we assume that the above analytical description of 
the evolution of E, 
pulse. Eq. 



p is also valid for the rise phase of the 
then gives 



t 

f 



1 



Eq 

E 



r,-l 



1 



(15) 



Here, Eq represents the maximal value of i?p, that is, the 
value at the start of the rise phase in contrast to the def- 
inition in Sect. 13.11 This equation is also valid for all 77, 
since as tends to 1, Eq. 1)15(1 has a limiting value of 



t 
T 



In- 



Eq 



(16) 



If the spectral evolution of a pulse is governed by this law, 
then the spectral lag. At, between two energy channels 
with characteristic energies Ei and E2, will be given by 



At 
~T 



Eq 
E2 



1 



E2 

El 



-n-i 



1 



(17) 



Here El < E2 and a positive lag is defined as the hard 
radiation precedes the soft. The ratio {E2/E1) in this 
equation is an instrumental constant and we denote the 
last two factors in Eq. IjlTfl by ■(/;(?? — l)j which tends to 
tjj{0) = \n{E2/Ei) as tends to 1. For a specific instru- 
ment the lag can thus be written as 



At^T-{EQ/E2r-' -^Kti-I). 



(18) 



Apart from the linear dependence on T, the lag depends 
on the value of Eq relative to the high energy channel, as 
well as to the instrumental function ■0(77 — 1). In Eq. ((18|l 
one could have chosen to include the second factor into 
the instrumental function. However, the energy band of a 
certain instrument limits the observed values of Eq, which 
therefore will be correlated with the typical energy of the 
band. Therefore, it is useful to keep this factor separate 
in such a comparison. Finally, using the definition of T, in 
Eq. ©, we find that 



^0 



ET 



which, as 77 tends to 1, has a limiting value 



At = $0 



Eq_ 

Eq 



V(0), 



(19) 



(20) 



In Fig. m the function -0 is plotted versus the HIC 
power-law index 77, for four different detectors, BATSE, 
Beppo-SAX, FREGATE and Swift. For BATSE, 7/> varies 
only weakly while for the other three instruments it 
varies somewhat more. For instance, as 77 changes from 
1.5 to 2.5, varies by a factor of 3.2 (BATSE), 5.3 
(FREGATE), 7.1 (Swift), and 19 (Beppo-SAX), respec- 
tively. The typical energies used were E1/E2 — 32.5/200 
(BATSE), 18.5/215 (FREGATE), 5.1/82.5 (Swift), and 
7.5/370 (Beppo-SAX). The variation of V for BATSE 
is relatively small, while spectral lags vary over several 




Fig. 1. Instrumental-dependent factor, ^, as a function of 
the HIC index 77, according to Eq. JHl- The solid line is 
for BATSE, the dotted line is for FREGATE, the dashed 
line is for Swift, and the dot-dashed line is for Beppo-SAX 



orders of magnitude l|Norris. Marani. fc Bonnell 1. hood . 

Furthermore, Ryde & Svensson (2002) showed that Eq 
cluster narrowly in < Eq >^ 332 ± 45 keV for BATSE 
bursts. The conclusion is, therefore, that the lag is de- 
termined mainly by the time-scale of the light curve. 
However, this relation is probably weaker for Beppo-SAX, 
FREGATE and Swift bursts, since they have a stronger 
dependence on 77, which will cause a larger dispersion. This 
is the case in particular for the large energy interval used 
in Beppo-SAX. 

An alternative, compact derivation of the above, con- 
taining the same information is the following: 



At ^ J X^dE = - j ^dE ^^Q j (kE")-^ dE 



$0 



k(l-v) 



El-" - E 



.ri-l 



(21) 



where we first used the HFC (Eq. |2]) and then used the 
HIC (Eq. n). 

3.3. Improvements of the Analytical Model 

In the discussion above, the change in the channel light- 
curves is assumed to be captured by the Ep decay, and the 
change in flux is assumed not to alter the lag significantly. 
The full description is, however, given by 

i^ch(t) = -Fboi(i) / B(E,t)dE 

Jch 

= E,,oi(t)Ep(t) j '^^'^ B(x)dx (22) 

Here, B(E, t) describes the instantaneous spect rum, for 
instan ce, by the empirical GRB function (Band et alJ 
l|l993(l 1 and Fboi(t) is the bolometric energy flux, that 
is, the flux integrated over the whole spectrum. As men- 
tioned above, we assume a monotonically decaying Ep and 
in the last step in the above equation we note that the en- 
ergy dependence of B(E,t(Ep)) is normally in the ratio 
X = E/Ep. Note that Bdx = 1/Ep. The full spectral 
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evolution is easily treated numerically by using Eq. H22I) . 
and this will be done in Sect. ^ These results are later 
compared with the above, simplified analytical model. 

There are two further issues that could be taken into 
account. First, in the above analytical treatment the low- 
energy spectral slope is assumed to be constant. This is 
indeed the case in many pulses, see for instance Figs. 3 
and 5 in iRvde fc SvenssonI l|1999(l . However, in other cases 
the slope can change significantly, often with a soften - 
ing trend. This can be seen Fig. 2 in lCrider et al.l (IfJOT*). 
Second, the pulse could actually be a tracking pulse which 
would need a more elaborate analytical approach. Finally, 
the Ep{t) behaviors used in the above derivations are the 
most commonly observed. However, it could be useful to 
study other behaviors. However, all these issues will lead 
mainly to minor corrections of the main behavior and they 
will be addressed elsewhere. 

In summary, the first conclusion that can be drawn 
from the analytical model is that the time lag measures 
a combination of all the parameters describing the spec- 
tral evolution, $0, ?7, Eq and Fq. Specifically, Ea. H18(l shows 
that the lag has a linear relation to the pulse decay time 
scale, T. In an alternative formulation there is also a lin- 
ear depe ndence on <I>o, as in Ea. ((T^ . iKocevski fc Lianj 
1 2003albl ) found that the l ag is correlated to the $0, see 
also discussion in lSchaefeil l|2004l) . A stronger linear corre- 
lation should emerge between At and ^q/Fq. If ■0(?7 — 1) 
and the ratio Eq / E2 varies only slightly, such a correlation 
will emerge. 

4. Numerical Simulation 

The above analytical treatment of the spectral evolution 
is useful, even though a few approximations had to be 
made. Here, we will numerically simulate burst pulses us- 
ing the full description in Eq. (|22(l and measure the spec- 
tral lags. We will use the knowledge gathered from the 
high-resolution spectroscopical investigations, described 
in Sect.O 

4.1. Numerical Setup 

Kocevski, Ryde, & Liang (2003) expanded the analytical 
study of the pulse shape, made above, to also include the 
rise phase. We review this treatment here as we will make 
use of their result in the numerical simulations. As the 
evolution of Ep is still from hard to soft during the main 
part of the rise phase of a pulse, the energy flux and the 
Ep are anti-correlated during this phase. By arguing from 
physical first principles lKocevski et al.l l)2003(l studied sev- 
eral analytical shapes for the whole pulse that includes the 
rise phase and asymptotically approaches Eq. (jl2|l in the 
decay phase (which is the behavior we are interested in 
here). In most physical models both the peak of the en- 
ergy spectrum and the luminosity are proportional to the 
random Lorentz factor of the shocked electrons to some 
power: Ep{t) cx T^{t) ■ g{t), and F^it) oc T\(t) ■ h(t). The 
functions g{t) and h{t) parameterize the unknown time 



dependence on particle densities, optical depth, magnetic 
field, kinematics, etc. The correlation between hardness 
and the energy flux can thus be described as 



F^{t)^El 



m 



(23) 



with 7 = b/a and f{t) = h{t)g{t)^^ . During the decay 
phase of a pulse the power-law relation dominates the 
HIC, as in Eq.^, whereas during the rise phase it will 
be dominated by the unknown function f{t). It is reason- 
able to assume that the rise phase is connected to some 
transient process, for instance, an initial decrease in opti- 
cal depth, increase in the number of energized particles, or 
the merging (or crossing) of two shells. After this initial 
phase, the original decay behavior as described by Eqs. 
(I12|l and H13|) should emerge. We therefore use the follow- 
ing prescription 



fit) a 1 - e-*/*^ 



(24) 



where the time constant ti- now represents the rise phase. 
This corresponds to h{t) cx g{t)^ a.s t >> t^ and the func- 
tion reaches the power law in Eq. asymptotically. A few 
new parameters have to be introduced; the energy flux, the 
peak energy, and the time of the peak of the light curve, 
^m, E^, and 



F = F^- 



Ep 
E„ 



1 



= -tm/ir 



(25) 



which, combined with the HFC (Eq.PI), gives the dif- 
ferential equation governing the spectral and temporal 
evolution of th e pulse . It has the following solution (see 
iKocevski et"al] l|2003|) for details). 



Fit)^ 



Ao(l 



-t/tr 



{l + {t + t,e-*/t')/T)'' 



{l + {t + t,e-*/t')/T)('^-^') 



(26) 



(27) 



where ^0 E^nd Ai are analytical functions of 
d,T,t^,^Q, F^, E^ and t^ and the decay index d is 
defined by requiring that F{t) t~'^ as t ^ 00, which 
gives r] = d/{d — 1). We introduce the peak energy at 
t^O.EQ which gives Ai = EQ{l + t,/T)^^/'^-'^\ The value 
of the flux at t = is arbitrary and is not import ant. The 
parameter d is found to be about 2-3 (,Kocevski et all 
2003D. 

The instantaneous photon spectrum is modelled by the 
standard GRB model ( Band et al. (, 1993,') ). which is essen- 
tially a low-energy power-law, F oc E'^ exponentially con- 
nected to a high-energy power-law F oc E^ where the pho- 
ton indices /3 < a. In an EF^ representation the spectrum 
peaks at Ep. The value of /3 was kept constant throughout 
the analysis here, since we concentrate our efforts to study 
the dependence of the lag on a. 

To be able to study spectral lags, the light curves for 
the four BATSE channels were found by integrating the 
spectra over the respective bands. Normally, the energy 
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Fig. 2. (a) Spectral lag as a function of pulse decay time scale, r, showing near-to linear relations. The parameters 
determining the spectral evolution were Eq = 500 keV, tr = I s, a = (black lines) and a = —0.5 (grey lines). The 
HIC power-law index (dlogF/dlogEp) rj = 1.3 (solid), rj = 1.5 (long dashed), 77 = 1.7 (short dashed), 77 = 2.1 (dotted), 
(b) Spectral lag as a function of rj. Here Eq = 500 keV, <r = 1 s, r = 5.5 s (black) and r = 2.5 s (grey). The low-energy 
power-law index, a = 0.5 (solid), a = 0.0 (long dashed), a = —0.5 (short dashed) and a = —1.0 (dotted). 
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Fig. 3. Contour plots of the spectral lag as a function of the HIC power-law index, r\ and the spectral, low-energy, 
power-law, a for varying values of Eq: (a) Eq = 100 keV, (b) Eq = 250 keV, (c) Eq = 500 keV, (d) £0 = 750 keV. 
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edges vary somewhat from observation to observation. For 
the simulation we used typical, average values: Ei — 24 
keV, E2 = 60 keV, E3 = 108 keV, E4 = 10000 keV. 
The spectral lag was then found by measuring the time 
between the peaks of the light curves in channel 3 and 
channel 1. The parameters that are used as input to the 
numerical analysis are the following: the maximal peak 
energy Eq, Fq, the low- and high-energy power laws, a 
and /3, the HIC index ry, the decay- and rise- time scales, 
T and t].. 

4.2. Spectral lags 

The solid lines in Fig. [21i show the spectral lag. At, as a 
function of the decay time scale t with a = 0, and for four 
different values of 77 = 1.3,1.5,1.7, and 2.1. The dashed 



lines show the corresponding behavior for a = —0.5. For 
both cases, Eq — 500 keV and ir = 1 s. The relation be- 
tween At and t is approximately linear (except for the 
case of 77 = 2.1, where a maximal value is reached). This 
shows that the lag is closely related to the decay timescale 
of the pulse. This is in agreement with the analytical re- 
sults in the previous section. 

The dependence of the spectral lag with the HIC 
power-law index, ry, is shown in Fig. [2h, for four differ- 
ent values of a = 0.5, 0.0, —0.5, —1.0. The solid curves are 
for r = 5.5 s and the dashed are for r = 2.5 s. Again 
Eg = 500 keV and t^ = 1 s. The lag has a maximal value 
for a certain 77- value, which increases for increasing a. This 
illustrates the complexity of the interpretation of the spec- 
tral lag: a certain lag does not uniquely correspond to a 
77 value. The longer decay time compared to the rise time 
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does not change the general behavior, it merely moves the 
curves to larger lags and larger 77- values. 

A more general picture of the dependence is given in 
the two following figures. In Fig.Olthc contour plots of At 
is shown as a function of the HIC index, 77, and a. The 
largest values are found for hard spectra (large a) and 
averaged HIC slopes. An increase in Eq, increases the lag 
values mainly for the hardest spectra. Fig. [Sja represents 
slices of these contours. The corresponding plots with a 
dependence on the maximal energy, Eq is shown in Fig. 01 
for a = 0.5,0.0, -0.5 and -1.0 (left to right). Both Figs. 
121 and 01 are for a pulse with = 1 s and t — 5.5. As a 
becomes softer (decreases) the lag, for a certain 77 and Eq, 
decreases. Also the maximal value occurs for smaller and 
smaller 77 values. We have made similar runs with varying 
pulse asymmetry, that is, varying the ratio of the rise and 
decay time-scales. As the decay time-scale, r, increases 
relative to t^-, the contour plots have in general a similar 
form except near the maximal value of the lag, that is, the 
peak of the contour plot. The peak increases and becomes 
more accentuated while the parameter space away from 
the peak is constant. 

4.3. Hardness Ratios 

We also studied the dependence of the hardness ratio (HR) 
on varying Eq, 77, and a over a single pulse. We simu- 
lated pulses and calculated the integrated photon flux in 
BATSE channels 1 and 3, over the interval of the pulse 
which was brighter than 0.1 of the peak flux. The ratio 
of these fluxes is defined as HR31. The results for a pulse 
with a rise time of 0.5 s and decay time of 2 s are depicted 
in the following three figures. Panel (a) in Fig.|31shows the 
contour plots of the HR as a function of 77 and Eq while 
a was fixed and set to 0.0. For low values of Eq the HR13 
is practically independent of 77 and is thus mainly deter- 
mined by Eq. This is the case in particular for large 77. A 
corresponding behavior is seen in panel (b) which shows 
the HR31 contours as a function of a and Eq. Here 77 was 
fixed instead and set to 2.0. The hardness ratio HR13 is 
practically independent of a for low values of Eq, most 
clearly pronounced for the hardest spectra. The strongest 
dependence of the HR is therefore on Eq. 

Often the relationship between the HR and the spec- 
tral lag is investigated and s uch plots are often used in 
LRS w ork (see, e.g.. Fig. 4 in lNorris. Marani. fc Bonnell I 
l|2000')). In panel (a) in Fig. the HR31 is plotted ver- 
sus lag for a run with a set to 0.0. The data points are 
connected for constant £^0, with Eq = 1000 for the curve 
furthest to the right and Eq — 100 for the curve furthest 
to the left. The index 77 increases upwards in the figure, 
along the curves, from 77 — 1.2 to 77 = 3.1. For low Eq the 
hardness ratio is seen to be practically constant, indepen- 
dent of 77, consistent with Fig. 5a. For large 77-values the 
dependence of HR on lag becomes greater (the vertical 
section of the curves). The large change of the lag as 77 
varies for a constant £^0 is also caught by Fig. 4b, as dis- 
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Fig. 7. The dependence of the hardness ratio HR'dl on 
the peak energy Eq, for the same run as in Fig.[S^. Each 
curve is for a different 77, increasing upwards in the figure, 
from 77 = 1.2 to 3.1. A variation in 77 will thus lead to a 
dispersion, which becomes larger for higher i?o-values. 



cussed above. In panel (b) the HR31-spectral lag relations 
are shown for varying a with 77 set to 2.0. The relations are 
approximate power-laws. This is consistent with previous 
figures; in Figs. 3, for a constant 77 = 2 and Eq, the lag has 
a monotonic increase as a function of a, while in Fig. 
again with j] — 2 and a constant Eq, we see that the HRs 
as a function of a will have a monotonic increase, or will 
be constant (independent) . Hence the correlations in panel 
Fig. Eb- Here, a larger Eq corresponds to a power-law at 
higher HR31 and the spectra get harder (larger a) to the 
right in the figure. Again, for low Eq values, the hardness 
ratio is independent of variations in lag. The HR31 and Eq 
are cor related which can be seen in Fig. ]^ and b. Among 
others. iBand et all 1 20041) have shown that such a relation 
exists observationally between Ep and the hardness ratio. 
This is shown in more detail in Fig. [7| which is from the 
same run as Fig.jS^; 77 increases upwards, from 77 = 1.2 to 
77 = 3.1. The dispersion, introduced by a variation in 77, is 
thus largest for large values of Eq . 



5. Discussion 

The above numerical work has presented the somewhat 
complex expression that the spectral evolution has in the 
two observables, namely the spectral lag and the hardness 
ratio. The results can be used to interpret lag and HR 
observations of various burst samples. 

5.1. Interpretation 

The dependence of the lag on a is given in Fig. 3. It shows 
that the lag increases with a for all tested values on 77. 
The more energetic the pulses/bursts are (higher Eq) the 
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Fig. 6. Hardness ratio as a function of spectral lag. In panel (a) each curve is for a certain Eq, with the thick line 
representing the largest value; rj increases upwards along the lines. Similarly, in panel (b), Eq increases upwards for 
each curve and a increases to the right along each curve. In panel (a) a = 0, and (b) r] = 2.0. 



steeper this correlation will be. Furthermore, the size of 
the lags can give a hint of the hardness of the spectra: a 
hard spectrum, for instance a thermal spectrum (a = +1) 
or a small pitch-angle synchrotron-spectrum (a = 0), will 
have a larger lag than a corresponding softer spectrum, 
maybe produced by optically-thin synchrotron emission 
{a = —2/3,-3/2). Figure 3 also shows that, if a sample 
of spectrally hard pulses (large a) is chosen, then the plot 
between lag and r] will have a characteristic bump indicat- 
ing the distribution of Eq in the sample studied. Note that 
each panel in Fig. 3 corresponds to a specific ii^o-value. 
A softer sample will have a less pronounced dependence 
(compare Fig. 2b). 

Furthermore, by studying the lag as a function of En, 
for a sample of bursts, conclusions can be drawn using 



the information available in Fig. 4. For instance, if the 
lag is found to increase with Eq, the sample is probably 
dominated by hard pulses, or at least do not have partic- 
ularly large ry-values. Oppositely, the lag could be found 
to decrease with Eq- This could indicate that the sample 
is dominated by soft spectra (panel d; especially at high 
energies). For harder spectra, such a decline is foimd only 
in pulses with the very largest rj. Continuing with Fig. 4, 
we see that if, once again, lag is plotted versus r] and only 
a weak dependence is found, the pulses are most probably 
low-energetic bursts, such as X-ray flashes. This is so es- 
pecially in the case of soft spectral pulses (compare Fig. 
3). For high-energetic bursts (high Eq) the characteristic 
dependence, with a peak at some intermediate rj, will be 
seen and the range of observed lags will be large. For very 
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soft pulses (small a) it will mainly be a declining func- 
tion. Note that the panels in the figure are for constant 
a-values. Correlations between, for instance, rj and a in 
the sample under investigation will complicate the inter- 
pretation slightly, since the same panel cannot be used for 
all rj. Also, since there are different possible dependence, 
the distribution of Eq of the sample will be important in 
determining the final relation. 

Similarly, differences in HIC-index rj between samples 
will be marked by that the hardness ratio generally is 
largest for the large-77 pulses (except for very low Eq ) . This 
is for instance manifested in Fig. 7 in which the large-r; 
pulses have a steeper dependence on the peak energy. The 
largest lags are found for pulses with intermediate 77 ~ 1.5 
- 77 ~ 2. As mentioned above, soft burst will have a char- 
acteristic switch in lag versus Eq behavior, from increasing 
(for low 7]) to decrease (for large 77). 

Figures 5, 6 and 7 reflect the behavior of the hardness 
ratios. If, for instance, the HR is plotted versus a (rj) and 
the relations are found to be weak, it is again an indication 
of low-energetic bursts. The distribution of the pulses in 
the HR-lag-plane will give indications of both Eq and rj 
according to the plot. 

5.2. Radiation processes and emission sites 

Several different emission mechanisms are probably active 
during the burst (see e.g. Ryde 2004), most importantly 
various versions of synchrotron (and/or inverse Compton) 
emission, predominantly its optically-thin version, and 
thermal, black-body emission. Other possible radiation 
mechanisms include synchrotron emission from electrons 
with a small pitch-angle distribution (or jitter radiation) 
and saturated Comptonization. The most visible differ- 
ence is in the hardness of the spect measured by the 
low-energy power-law index a. A certain sample of bursts 
can thus be studied with the relations discussed above, to 
discern the dominating radiation process. For instance, on 
the basis of the durations and hard nesses, burs ts seem to 
belong to two distinct classes e.g. iBalazs et all ll2003l|) '). 
or even three l|HorvathL Il998l: iHorvath et all 12004 - Anv 
difference in radiation mechanism between these classes 
can therefore be efficiently investigated. 

It is also known that the importan ce of the curvature 
effect varies between pulses (see e.g. iRvde fc PetrosianI 
In the case of it being important, the HIC index, 
77, will have a characteristic behavior and a large value. In 
Fig. [SJl, 77 = 3 for a curvature-driven pulse is marked by 
the dashed line and the relation between IIR31 and spec- 
tral lags for such pulses is shown in Fig.|Hl This shows the 
hardness ratio decreasing linearly with the corresponding 
time lag. In this particular run a — —0.5. Eq decreases to 
the right with the point furthest to the left with Eq ~ 1000 
keV and the point furthest to the right with Eq = 100 
keV. This decrease is obvious in Fig. 4, in which large val- 
ues of 77 and soft spectra sample a space where the lag 
decreases with Eq. This, together with the positive cor- 
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Fig. 8. Correlation between hardness ratio and spectral 
lag for a "curvature-driven" pulse (77 = 3). Eq decreased 
to the right. See the text for further details. 

relation between Eq and HR31 depicted in Fig. then 
gives the relation in Fig. |H1 On the other hand, if the 
dominant time-scale is not the curvature time-scale then 
the intrinsic 77, representing the comoving radiation pro- 
cess, will be revealed. For instance, assuming synchtroton 
emission and a constant bulk Lorentz factor, F, we expect 
Ep oc 7,^„j„ • B (the minimum electron Lorentz factor and 
the magnetic field strength) and F (x £e ■ ■ ^min, where 
£e ~ /q°° jN{j)d'-f is the total electron energy (^^(7) is 
the energy spectrum of electrons). If the latter quantity 
is constant, we expect rj — 2 ii B varies and 7min is con- 
stant and 77 = 1/2 if i? is constant and 7min varies. The 
observed correlations would then sample different regions 
in the panels discussed above. 

During the past few years increasing numbers of X-ray 
rich GRBs and so-called X-ray flashes (XRF) have been 
observed (Heise et al. 2000, Barraud et al. 2003). Several 
studies have shown that the main difference in properties 
between these populations and classical GRBs is the peak 
energy distribution. It has therefore been suggested that 
they are all of the same origin (e.g. Kippen et al., 2002, 
2003, Lamb et al., 2004) and that XRF is an extension 
of the GRB energy distribution. As described above the 
distribution of peak energy between different samples will 
have characteristic signatures on the lag correlations and 
can be further used in studying the difference between 
GRBs and XRFs. 



6. Conclusion 

The two spectroscopical analysis methods discussed in this 
paper, HRS and LRS measure the same underlying be- 
havior and are alternative ways to represent the spectral 
evolution. We have shown analytically and numerically 
that the spectral evolution described by high spectral- 
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resolution data, in empirical correlations, naturally leads 
to correlations involving spectral lags. The interpretation 
of the spectral lag is not straightforward. This depends 
on the fact that the lag measures a combination of spec- 
tral parameters such as r], Eq, and r (apart from effects 
caused by a and /3; their values and evolution). It fur- 
ther depends on the relative channel widths into which 
the data are divided. For different pulses within a burst, 
the spectral parameters are in general not the same, and 
thus the lag is not expected to be the same during the 
whole burst. This is indeed the case in observed bursts, as 
we note in paper II (Ryde et al. 2004). 

The main result of the above analysis is that the lag 
correlates strongly with the decay time-scale of the pulse 
(or equivalently, of the peak-energy decay) and that this 
relation is linear. This means that the relationships includ- 
ing the lag, that have been identified, should be translat- 
able into relationships including the pulse time-scale. This 
time-scale is closely connected to the processes forming 
the pulse in the out-flowing plasma; the radiation pro- 
cesses, acceleration mechanisms, light-travel effects and 
outflow velocity etc. In particular, we find a good correla- 
tion between l ag and the quantity F/ $n , reproducing the 
observation bv lKocevski fc Liangl l|2n08bj) . We further find 
that harder spectra, with large a, have the largest spectral 
lags. This is, in particular, the case for averaged 77-values. 
Similarly, an increase in Eq, increases the lag, except for 
soft a (see Fig. 4 d). Interestingly, the hardness ratio is 
found to be only weakly dependent on rj and a, and in 
particular for low, initial peak energies Eq it is practically 
independent, and it is mainly determined by Eq; according 
to Fig. 6a. The relation between the hardness ratio and the 
lags, for a certain Eq and ry, are described by power-laws, 
as a varies (Fig. 5d). Furthermore, for a curvature-driven 
pulse {rj = 3), the hardness ratio decreases linearly with 
lag for lower and lower Eq- 

The analytical and numerical results presented here 
are the consequences of the empirical description of the 
spectral evolution in pulses, which have been firmly es- 
tablished. These results can thus be used as a reference 
in the analysis of observed pulses and bursts, in partic- 
ular low spectral-resolution studies in which the spectral 
parameters cannot be found and mainly the lag and hard- 
ness ratios are analyzed. Such an analysis is performed 
in paper II (Ryde et al. 2004) in which we also discuss 
how the conclusions presented here relate to the physical 
models that have been proposed to explain the observed 
lag-correlations . 
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